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SUMMARY 

An analysis and discussion is presented of existing estimates of 
the J-integral for cracks in infinite bodies. Equations are presented 
which provide convenient estimates for Ramberg-Osgood type elasto- 
plastic materials containing cracks and subjected to multi axial 
loading. The relationship between J and the strain normal to the crack 
is noted to be only weakly dependent on state of stress. But the 
relationship between J and the stress normal to the crack is strongly 
dependent on state of stress. A plastic zone correction term often 
employed is found to be arbitrary, and its magnitude is seldom 
significant. 


^Material also submitted to Engineering Fracture Mechanics for 
publ ication. 
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INTRODUCTION 


Materials are considered which have a uniaxial stress-strain 
relationship of the form: 

e/e 0 = o/o 0 + a(o/o Q ) n (1) 

where o is stress, a , e , and a are constants more fully explained 

o o ■ 

later, and n is a constant generally called the strain hardening 
exponent. The first and second terms represent, respectively, elastic 
and plastic components of the total strain, e. 

Configurations of cracks in infinite bodies which are of interest 
are shown in Fig. 1. Two dimensional cases as in Fig. la may be either 
plane stress or plane strain, depending on the situation in direction 
3. Note that a stress, S, is applied remote from the crack in a 
direction normal to the plane of the crack. Remotely applied transverse 
stresses, T, are also considered. Figure lb illustrates remotely 
applied axisymmetric loading of a circular crack. The stress normal to 
the plane of the crack is again S, and T is the stress in all radial 
directions in the plane of the crack. Edge or surface crack 
configurations correspond to bisecting either Fig. la or lb along the 1- 
3 plane. 

Some early. work related to this discussion is that of Shih [1] and 
Shih and Hutchinson [2], which established the general form of estimates 
of J for cracked members of Eq. 1 type material. The EPRI Elastic- 
Plastic Fracture Analysis handbook [3] later pursued such estimates in 
considerable detail, but only for cracks in finite bodies. He and 
Hutchinson [4,5] have reported analytical results for infinite body 
cases as in Fig. 1 that figure prominantly in what follows. They 



considered multiaxial loading (various ratios T/S) and remote strains 
that are predominantly plastic, so that Eq. 1 reduces to: 

e/e o = a(o/o 0 ) n (2) 

Numerical results for certain cases of uniaxial loading (T = 0) are 
also reported by Trantina et. al . [6,7] using Eq. 1. In addition, 

Wilson [8] has previously studied two-dimensional cases of uniaxial 
loading, considering not only Eq. 1 but also bilinear and piecewise 
power hardening. Also, I previously presented limited discussion of 
this subject in Appendices to two papers [9,10] that included use of 
such J-integral estimates in analyzing growth data for small fatigue 
cracks. 

It is the purpose of this paper to provide a useful analysis and 
discussion of J-integral estimates for cracks in infinite bodies of Eq. 

1 type materials. The estimates are reduced to a convenient form for 
application to engineering or research problems. Also, discussion is 
presented concerning certain important aspects of the estimates, 
including multiaxial effects, and a plastic zone correction term. 



STRESS-STRAIN RELATIONSHIPS 

Before proceeding, it is useful to further explain the uniaxial 
stress-strain relationship, Eq. 1, and also to extend this to cases of 
multiaxial loading. 

Uniaxial Stress-Strain Curve 

Figure 2 illustrates a curve corresponding to Eq. 1. The curve 
gradually diverges from an elastic line having a slope: 

E - o 0 A 0 (3) 

Hence, (e Q , a ) is a point on the elastic line, not on the actual 
curve. The constant o Q may be thought of as a yield strength defined in 
the conventional manner by a plastic strain offset. Then (e , o Q ) is a 
point on the actual curve, where 

= (1 + o)e 0 (4) 

And the plastic strain offset is: 

e = a t 

py o 

Of the three constants e , o , and a, any one of them may be 

arbitrarily chosen, and then the other two fitted to the particular 

curve of interest. For a given set of values of these constants, there 

is a family of curves, all passing through the point (e , o Q ), and each 

corresponding to a different value of n. Hence, the value of n also 

needs to be fitted to the particular curve of interest. 

Of the three constants e , o , and a, the latter is the most 

o’ o 

logical to be arbitrarily chosen. The constant a can be thought of as 
representing the relative magnitude of the plastic strain offset. 
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a 



0 



( 5 ) 


Conventional plastic strain offsets, such as e p ^ = 0.002 for engineering 
metals, generally correspond to a in the range 0.1 to 1.0. 


Multiaxial Stress-Strain Relationships 

All of the equations to be used later involve only the stresses and 
strains remote from the crack. No remote shear stresses are applied and 
the material is assumed to be isotropic, so that the axes shown in Fig. 

1 are also the principal axes for both stress and strain. These 
principal stresses and strains will be denoted by a. and e., where i 
corresponds to any one of the directions in Fig. 1. Subscripts j and k 
then indicate the remaining two directions. 

Strains may be separated into elastic and plastic components. 

i ei pi v ' 

where the subscripts e and p indicate these components. The stresses 
and the elastic strains are related by Hooke's law. 

e . = I- [o. - v (o . + o, ) | (7) 

ei E i v j k' v ' 

where v is the elastic constant, Poisson's ratio. (A value of 
v = 0.3 is used throughout this paper as being typical for engineering 
metals) . 

As only constant ratios T/S are considered, proportional remote 
loading applies, and a total deformation plasticity theory can be used 
for the remote stresses and strains. In this context, stresses and 
plastic strains are related by: 

e p . = [o i - 0.5 (oj + a k )] (8) 
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where 


e p = 5/e p 


(9) 


The quantity o is the effective stress 


r (°i " ° 2 ) 2 + (°2 " °3) 2 + (°3 " CI l) 

0 = l 5 J 


2 1/2 


( 10 ) 


and ep is the effective plastic strain 


" = I — (c 2 + e 2 + e I 1/^ 

p l 3 { p> p 2 p’ ;| 


(11) 


Substituting Eqs. 7 and 8 into Eq. 6 gives the following 
relationship between stresses and total strains: 


E . = 




- V 


(°i + ] 


where 


E. = a/e 


(12) 


(13) 


The quantity e is the effective total strain, given by: 


e = a/E + e 

P 

and v is an effective Poisson's ratio. 

va + 0.5 Ee 

y = — 

Ee 

Note that as s is increased from predominantly elastic to predominantly 
plastic strain, v increase smoothly from v to 0.5, so that Eq. 12 
reduces to Eq. 7 in the former case, and to Eq. 8 in the latter case. 
Also, substituting the inverse of the three equations represented by Eq. 
12 into Eq. 10 gives an expression for e as a function of the principal 
total strains: 


(14) 


(15) 
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(e t - e 2 ) 2 + ( 


e 0 - e 


i ) 2 + (e 3 - e t ) : 


1/2 


e = 




(16) 


.1 + 5 

For predominately plastic strains, where v = 0.5, Eq. 16 reduces to Eq. 
11. And for predominately elastic strain, where v = v, substitution of 
Eq. 7 into Eq. 16 gives Eq. 10. 

Effective stresses and plastic and total strains as defined above 
coincide with the corresponding quantities for the uniaxial case, so 
that Eq. 1 can be generalised to: 


e/e Q = o/a Q + a (a/a Q ) n (17) 

Hence, the constants e Q , o , a, and n for use in Eq. 17 can be evaluated 
from a uniaxial stress-strain curve. 

For the particular cases of interest of Fig. 1, Eqs. 10 and 12 may 
be used to derive ratios between the effective stresses and strains and 
the corresponding stresses and strains in the direction normal to the 
plane of the crack. These ratios, which will be needed later, are given 
below. The loading in the plane of the crack is characterized in terms 
of the quantity R. 

R = T/S (18) 

(a) Plane stress: Oj = S, a 2 = T, a 3 = 0 

a/S = /I - R + R 2 (19a) 

e/e 1 = — (19b) 

1 - vR 


(b) Plane strain: o t = S, o 2 = T, e 3 = 0 


a/S = /(I + R) 2 (1 - v + 5 2 ) - 3R 


(20a) 
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(c) Ax i symmetric case: 



RESULTS OF HE AND HUTCHINSON FOR ELASTIC STRAINS SMALL 

In the work of He and Hutchinson (4,5) numerical results for the J- 
integral, Jp, for Eq. 2 type material are presented by tabulating values 
of hj = h, (n,R) , 
where 

J 

h t (n,R)=^— (22) 

oea 

Various combinations of n and R = T/S were investigated for the 
configurations of Fig. 1. 

In addition, a perturbation approach was used |4,5| to obtain 
closed form estimates for the plane strain and axisymmetric cases. 

These are as follows, with the obvious guess for plane stress also being 
included: 

(a) Plane stress: 

J f.S 2 

_JL_ = h, 1 = it/n (— ) (23) 

oe:a a 

(b) Plane strain: 

J _ - f,S ' 

= h,' = ^ (— ) (24) 

oe a o 

(c) Axisymmetric: 

J c f.S 2 

-j— =h/= - 6 (— ) (25) 

aea n/1 + 3/n a 

The factor f t = l for central or interior cracks as in Fig. 1. It 
is introduced so that cases of edge or surface cracks are included, with 
these being described by bisecting Fig. la or b along the 1-3 plane. 
Values of f t will be estimated from the linear elastic case (n = 1), for 
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example, f, = 1.12 for two dimensional edge cracks formed by bisecting 
Fig. la. 

Figures 3-6 compare each numerical result from refs. 4 and 5 with 
the closed form estimate from one of Eqs. 23-25. Ratios of the 
estimated value h , ' to the numerical result h t are plotted for each 
value of the strain hardening exponent, n, where there is a numerical 
result. The numerical h t values are either upper or lower bounds, with 
both bounds being available in only a few cases. Arrows are therefore 
shown in Figs. 3-6 to indicate that the true ratio h^/h, may be above 
or below the value plotted, as appropriate in each case. 

Figure 3 gives the comparisons for central cracks under plane 
strain, where there are numerical results for T/S = 0, -1, and 0.5. 
Figure 4 gives comparisons for edge cracks under plane strain, and for 
central and edge cracks under plane stress, for which cases numerical 
results are available for T/S = 0 only. The ratio h^/h for these cases 
deviates increasingly from unity as n is increased up to the highest 
value used of 7. Most of the values are within 10% of unity, and all 
are within 16%. 

Note that both bounds are available for central and edge cracks 
under plane strain with T/S = 0. Here, h^/h, for the middle of the 
bounds is below but within 5% of unity for central cracks, and above but 
within 11% of unity for edge cracks. As pointed out in Ref. 5, this 
implies that f t is actually not constant at the value of 1.12, but 
decreases with n, perhaps approaching unity for large n. 

Figures 5 and 6 give h t ' /h 1 ratios for the wide range of T/S values 
studied for the axisymmetric case. Lower bounds only are available from 
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Ref. 4 for n up to 10, and also for n = . The figures show all but 

the n = ® results, which are as follows: 

T/S 0 -1.-4 -9 .5 .67 .75 


hj'/hj 1.037 1.073 1.087 1.093 1.146 0.942 0.549 

For T/S < 0, the ratio h 1 1 /h t is within 10% of unity for all n values. 
But for T/S > 0, the discrepancy is almost 50% for the combinations of 
large n and high T/S value. 

If it is assumed that the numerical bounds on h, for the various 
cases of Figs. 3-6 are reasonably close to the actual values, say within 
10%, then the simple estimates from Eqs. 23-25 are suitable for most 
purposes. Note that the Only large discrepancies occur for axisymmetric 
cases approaching triaxial tension (T/S = 1), which state of stress is 
rare in practical applications. Caution is obviously needed for large n 
values, and additional numerical results are needed to confirm or modify 
this tentative conclusion. 
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ESTIMATES FOR ELASTO-PLASTIC MATERIALS 

Based on the work of Shih (1) and Shi h and Hutchinson (2), J- 
integral values for Eq. 1 type materials may be estimated as follows: 

J = J e + J z + J p (26) 


where Jp is the estimate of J for Eq. 2 type material which has a stress 

versus strain curve that is the same as the stress versus plastic strain 

curve for the Eq. 1 type material of interest. (This would be satisfied 

if the constants e , o , a, and n are the same.) 
o o 

The term J g is calculated for a linear-elastic material with the 
same elastic modulus, E = ° 0 A 0 > as the Eq. 1 type material. 


V7 C nS 2 a 
i _ K * __ o 

J e = V = E 1 


(27) 


where 


E' = E 


(plane stress) 


E 1 = E/(l - v 2 ) 


(plane strain, and 
axi symmetric) 


The quantity K is the Mode I stress intensity factor of linear 
elastic fracture mechanics, which is given by: 

< = fjf 2 S/7a (28) 

Note that the geometry factor is separated into the product of a free 
surface correction factor, f t , as in Eqs. 23-25, and a residual geometry 
factor, f 2 , with values as follows: 
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Case 

f. 

f 2 

C o 

Fig. la - Central crack, plane stress 
or plane strain 

1.0 

1.0 

1.0 

Bisected Fig. la - Edge crack, plane 
stress or plane strain 

1.12 

1.0 

1.25 

Fig. lb - Internal circular crack 

1.0 

2/n 

0.405 

Bisected Fig. lb - Half-circular 

1.035 

2h 

0.434 


surface crack 

From comparison of Eqs. 27 and 28, note that: 

C Q = V V (29) 

The values of f t and f 2 are the well known values available from 
any good collection of linear-elastic fracture mechanics results, such 
as Ref. 11. As shown in Ref. 6 and previously elsewhere, f L for a half- 
circular surface crack varies around the periphery of the crack, 
decreasing from about 1.2 at the intersection with the free surface to 
about 1.035 at the point of maximum depth. 

Note that J cannot be separated into elastic and plastic terms, J g 
and Jp, in a formal mathematical sense. This summation as in Eq. 26 is 
only a strategy adopted to obtain an estimate. It provides correct 
values at the limiting case of linear-elastic loading, that is, 
negligible plastic strain, where J = J g , and also at the limiting case 
of large plastic strain, that is, negligible elastic strain, where J = 
Jp. 

The term J z of Eq. 26 is based on a plastic zone calculation and is 
used in Refs. 1 and 2 to smooth the transition between the limiting 
cases. Following Ref. 3, the linear-elastic solution, Fq. 27, is 
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modified by replacing the crack dimension "a" with an effective value, 


a e* 


where 


Co - 


$ = 



a = a + 
e 

<ir 

y 

(30) 

= J,_ 
C 2 TT 

■ ( n_1 ) r K ) 

l n+l ; l o Q J 


(31a) 

= 2 


(plane strain) 

(31b) 

= 6 

1 

(plane strain, and 
ax i symmetric) 

(31c) 
( 3 Id ) 

1 + 

(P/PJ 2 



The quantity P/P 0 is the ratio of the applied load to a limit or 

reference load based on the stress a . For cases of cracks in infinite 

o 

bodies, and considering multiaxial loading, the following substitution 
is appropriate: 


P/P = a/a 

O 0 


(32) 


The modified value of J e is: 


a + <tr 

J ' = J ^ 

e e a 


(33) 


and the increase of J ' over J is defined as J . 

e e z 

<i>r 

J = J 1 - J — — J 
z e e a e 


(34) 


Combining Eqs. 28, 29, 31 and 34 gives: 

J C , e 

T 1 - r 2 (£f ) 0 


( CT / <T ) 
v 0 


J Cj n+1" - i [ - f \ 

e o 1 + (a/'i ) 


(35) 
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It is convenient to normalize the estimated J from Eq. 26 to the 


elastic term. 


J/J e = 1 + J z /J e + J p /J e (36) 

Only Jp/J e remains to be evaluated. This is obtained by noting that the 
effective plastic strain is given by the second term of Eq. 17. 


ip/e o = «(5/o 0 ) n 


(37) 


Also, e in Eq. 22 must be interpreted as equal to this e p , as J p is 
estimated for an Eq. 2 type material with negligible elastic strain. 
Then combining Eqs. 22, 27, 36, and 37 leads to: 


= 1 _ _L 

J 1 J C ttE 
e e o 


h,a E 1 - 2 - n-1 

(I) (M 


(38) 


Recall that o/S is given by the appropriate one of Eqs. 19-21. The 
combination of Eqs. 35 and 38, when reduced to the limiting case of 
uniaxial loading (T/S =0), is similar to Eq. 18 of the paper by Wilson 
( 8 ). 

A simplified version of Eq. 38 which employs the estimates 
of hj from Eqs. 23-25 is as follows: 

J J - 

J “ 1 + _ 
e e 

where, consistent with Eq. 14: 


i +1 ^h 

J o E - 

e 


(39) 


E - n/E 
e 


(40) 


And the expressions for h Q are: 

h = /n 
o 


(plane stress) (41a) 
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3/n 

4 


(plane strain) 


(41b) 


h = 
o 


2/1 + 3/n 


(ax i symmetric) (41c) 


In addition, it will be shown below that J z /J e is often small, so that 
reasonable estimates are given by: 


t = 1 + h o r ^ < 42 > 

e E e 

It is sometimes convenient to employ an equivalent K value which is 
modified from the elastic value based on the J-integral. Using Eqs. 27 
and 42 gives: 

1/2 

K =/3F= K (lth 0 ^] (43) 

£ e 

Hence, from Eqs. 42 and 43, a correction factor can be applied to a 
J e or K value for the linear-elastic case, with this factor depending on 
the strain hardening exponent through h Q and the relative degree of 
plasticity, as expressed by the ratio of the effective plastic to 
effective elastic strain. 
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DISCUSSION 


Discussion follows under three subheadings. First, there is a 
discussion of the trends in J-integral values with state of stress. 

Next, there are comparisons with existing analytical results for elasto- 
plastic materials. Finally, the plastic zone correction term is 
discussed. 

State of Stress Trends 

A question of some importance is whether the J-integral has the 
potential of predicting fracture and fatigue of real materials where 
multiaxial loading affects the behavior (12). First, note that the 
elastic term, Eq. 27, is not affected by the in-plane stress, T, as it 
depends only on S. The plastic zone term, Eq. 35, is affected by T, but 
this is unimportant, as this term is seldom a significant fraction of 
the total J value. But the plastic term, such as the last term of Eq. 

38 or 39, is strongly affected by T. This occurs as a result of the 
effective stress and strain quantities which appear. 

Figures 7 and 8 illustrate the trends for the ax i symmetric case for 
various ratios of the stresses T/S. Figure 7 plots J versus the 
principal strain in the direction normal to the crack, and Fig. 8 plots 
J versus the principal stress, a t = S, in this same direction. These 
curves are based on typical values of n = 7 and \> = 0.3, a value of 
a = 1 is arbitrarily chosen, and the plots apply for any choice 
of o Q and e . The approximate form represented by Eqs. 39 and 41c is 
used, except that Eq. 38 with an interpolated numerical h, value is used 
at T/S = 0.75. This is done as the approximation is less than about 5% 
except for T/S = 0.75, where it would have been about 25'/.. Note that 
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Eqs. 21a and b were also employed in the calculations to obtain these 
curves. 

It is interesting that J is only weakly affected by state of stress 
when comparisons are made based on principal strain, but drastically 
affected for comparisons based on principal stress. In addition, the 
trend of higher J values for lower T/S values is at least qualitatively 
consistent with the empirical trends observed in fracture and fatigue 
response of materials. 

Consider a given material and a given choice among plane stress, 
plane strain, or the axisymmetric case. Equation 38 then gives a nearly 
unique relationship between J and either effective stress, o, or 
effective strain, e, which is independent of the stress ratio T/S. 
Furthermore, if the approximations of Eq. 41 are used, and if the 
plastic zone term is omitted by the use of Eq. 42, the relationships 
between J and e, and between J and o, are exactly unique. 

Comparisons With Elasto-Plastic Analysis of Trantina 

References 6 and 7 present J-integral results from finite element 
analysis using an elasto-plastic stress-strain curve of the form of Eq. 
1. The approximation involved in summing separately computed elastic 
and plastic terms to obtain J is not necessary. Hence, these results 
provide a valuable comparison with the estimates presented above. 

All of the analysis in Refs. 6 and 7 is for uniaxial loading. 
Reference 6 considers the important three-dimensional case of a half- 
circular surface crack, and Ref. 7 considers circumferential cracks 
growing from spherical and other voids. Figures 9 and 10 compare these 
analytical results with estimates from Eqs. 35 and 39. The various 
constants used in the estimates are as follows: 
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Fig. 9 estimates based on: 
circular crack 2D crack 


Fig. 10 estimates 
(circular crack) 


C 0 0.434 

C 2 6 

E'/E l/(l-v2) 

h Q ‘ 3/2/1 + 3/n 


0.434 

6 

1/C1-V2) 

3/n/4 


0.405 

6 

l/(l-v2) 

3/2/1 + 3/n 


For a half-circular surface crack, estimates made based on h Q for a 
circular crack give curve 1 in Fig. 9, which falls considerably below 
the analytical results. Note that this case is a complex three 
dimensional one, and the special situation existing for a circular crack 
under axisymmetric loading may not exist even approximately. The 
possibility exists that the behavior is closer to that for a two- 
dimensional crack. Hence, an alternate estimate is made which is the 
same except that h Q is taken from Eq. 41b, which applies for two- 
dimensional cracks in plane strain. This yields curve 2 in Fig. 9, 
which is in excellent agreement with the analysis. 

Note that the comparisons with the analysis of Ref. 6 depend on the 
position around the periphery of the crack. The above comparisons are 
all for the point of maximum depth, and the appropriate value 
of f t = 1.035 is used in normalizing the values from Ref. 6 to obtain 
J/J e for plotting in Fig. 9. 

Reference 6 also contains limited analysis for n = 20. Similar 
comparisons were made, and similar results were obtained. In 
particular, h Q for circular cracks (Eq. 41c) gives estimates which are 
low by more than a factor of 2 at e/e Q = 1.35. But use of h Q for two- 
dimensional plane strain (Eq. 41b) gives an estimate which is high by 
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only 8 % at z/t Q = 1.35. 

Figure 10 gives comparisons of analysis and estimates for a 
circumferential crack at a spherical void. The estimates are based on 
the assumption that this geometry is equivalent to a circular crack of 
radius a = R + a, where R and a are defined in Fig. 10. In other words, 
the crack tip is assumed to be sufficiently remote that there is no 
effect of the local stress concentration of the spherical void. 

Figure 10 shows reasonable agreement between estimates and analysis 
for n = 5, 10 and 20. The J values from analysis at the higher strains 
are 10 to 20% higher than the estimates. Note that, for linear 
elasticity (n = 1), the effect at a/R = 0.5 of the stress concentration 
of the spherical void is to elevate the value of f 2 by a factor |13| of 
about 1.04, where the interpretation a = R + a is made for Eq. 27, so 
that: 

C Q = (1.04 x 2 /tt) 2 - 0.438 

Since this new value of C Q is used in normalizing the analytical results 
of Ref. 7 to obtain J/J e , the effect of the stress concentration has 
been removed from the comparisons to this extent. The fact that the 
analytical values are still higher than the estimates may indicate that 
the effect of the stress concentration on J/J 0 is greater for the higher 
n values. 

In preparing Fig. 10, it was necessary to know the e Q value used in 

Ref. 7. This value is e = 0.005776 and was obtained from information 

o 

kindly provided to the author by Trantina. 

Note that the approximation using a simple circular crack of radius 
a = R + a is the long crack limiting case for the actual problem of a 
circumferential crack at a spherical void. For linear-elasticity, the 
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difference between this case and the actual solution becomes 
significant, say > 20% in K, around i/R = 0.08. For small 5 ,/R, the 
actual solution deviates drastically from the long crack limiting case 
and becomes dominated by the local stress field of the spherical void 
[13,141. Similar behavior is expected for n > 1, so that estimates 
similar to those just described could not be reasonably applied to the 
analysis of Ref. 7 for a/R = 0.1. 

To summarize, the comparisons of Figs. 9 and 10 provide 
encouragement that the proposed estimates are reasonable. However, no 
analytical results for an Eq. 1 type material appear to be available 
which allow the inaccuracies of the estimates to be isolated. The Fig. 

9 comparison is complicated by the three-dimensional nature of the half- 
circular surface crack, and the Fig. 10 comparison is complicated by the 
spherical void. 

Plastic Zone Correction Term 

In the estimates of J for an Eq. 1 type material provided by Eq. 38 
or 39, the elastic (first) term dominates at low strain, and the plastic 
(third) term dominates at high strain. The plastic zone (second) term 
makes its greatest relative contribution at intermediate values of 
strain. This is illustrated by Fig. 11, where J 2 /J computed from Eqs. 

35 and 39 is plotted versus strain for an interior crack under plane 
strain with T/S =0. A typical value of n = 7 is chosen, and four 
different a values are employed. The strains plotted are normalized to 
a common value of yield strain rather than the individual e Q values 
corresponding to each a. This common yield strain was chosen as 
the Ey value from Eq. 4 corresponding to a = 1, called c 1 . Strains 
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normalized to this value are calculated from strains normalized to e by 

o J 

noting that: 

e y ' = e Q (l + a 1 ) (a7a) 1/( "' 1) 

where a'=l is the chosen standard value. 

Two significant facts are observed from Fig. 11. First, the 

plastic zone correction term has the undesirable property of depending 

on the arbitrary coice of a. And second, for any reasonable choice 

of a, typically 0.1 to 1, the J z term is never a significant fraction of 

the total J. For the case of Fig. 11, the J z term is always less than 

10% of J, even for the extreme choice of a = 0.01. 

The dependence on a arises from the fact that, for a given stress- 

strain curve, different choices of a result in different values of o . 

o 

The size of the estimated plastic zone decreases if a Q is increased, 
which corresponds to increased a, with the result that J z /J e 
decreases. This can be seen by studying Eqs. 3-5 and 35. 

Plots similar to Fig. 11 were made for a wide variety of cases, and 
the peak J z /J values were noted for each case. The approximate 
estimates employing h Q in Eq. 39 were used, except that Eq. 38 with the 
numerical h ( value was used where this was known to gave a difference in 
the plastic (third) term exceeding 15%. Figure 12 summarizes the peak 
effect of the plastic zone correction as a function of n for uniaxial 
loading for both a = 0.1 and 1. The largest effect is for plane stress, 
and this never exceeds 15% for n up to 10. 

Figure 13 shows the effect of various degrees of in-plane loading 
on the axisymmetric case for n values up to 10. The plastic zone term 
is most significant for states of stress approching hydrostatic tension. 



approaching 30% for T/S = 0.75 and n = 10. For plane stress and plane 
strain with n up to 10 and T/S up to 0.75, the estimated value of J z /J 
seldom exceeds 10%. The exceptions are plane stress at T/S = 0, as 
shown in Fig. 12, and plane strain at T/S = 0.75, where J z /J e arounc * 
20 %. 

Accurate numerical analysis for Eq. 1 type material would be needed 
for critical cases where J 2 /J is large to determine if the J z /J 
correction term is necessary and to guide in the choice of u. However, 
at present it appears that the term is generally unnecessary, in 
addition to being arbitrary due to the effect of u. 
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CONCLUDING REMARKS 


Estimates of J from the relatively simple Eq. 42 appear to be 
suitable for most practical purposes. More specific estimates (Eq. 38) 
may be needed for axisymmetric loading which approaches triaxial 
tension. For nonaxi symmetric cases, such as two-dimensional plane 
stress and plane strain, caution is needed for high tension loads in the 
plane of the crack, and for large values of the strain hardening 
exponent, n, where analytical support for Eq. 42 is weak. 

The dependence of the plastic term of the J estimates on the ratio 
of effective plastic strain to effective elastic strain is 
significant. This corresponds to a dependence of J on state of stress 
which is quite large for a given value of stress normal to the crack. 

Additional work is needed on the complex three-dimensional, but 
practically important, case of a half-circular surface crack. 

Preliminary indications are that the J values exceed estimates based on 
the axisymmetric case. 

The plastic zone correction term is both small and arbitrary. It 
should be neglected until additional analysis establishes its need and 
removes its arbitrariness. 
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Fig. 1 - Cracked body configurations considered: (a) Two dimensional 

cases, either plane stress (a 3 = 0) or plane 
strain (e 3 = 0). (b) Circular crack under axisymmetric 

loading. Bisecting (a) along the 1-3 plane describes cases of 
two dimensional edge cracks, and bisecting (b) along the 1-3 
plane gives a half-circular surface crack. 
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cr, Stress, MPa 



Fig. 2 - Illustration of a uniaxial stress-strain curve. 
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- Comparison of J-integral constants from He and Hutchinson for 
central cracks under plane strain. The numerical hi values are 
either upper or lower bounds - see arrows. 
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Fig. 4 - Comparison of J-integral constants from He and Hutchinson for 
edge cracks under plane strain, and for central and edge cracks 
under plane stress. The numerical hi values are either upper 
or lower bounds - see arrows.. 
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Fig. 5 - Comparison of J-integral constants from He and Hutchinson for 
circular cracks under axi symmetric loading, either uniaxial or 
with in-plane compression. The numerical hi values are all 
lower bounds - see arrows. 
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Fig. 6 - Comparison of J-integral constants from He and Hutchinson for 
circular cracks under axisymmetric loading with in-plane 
tension. The numerical hj values are all lower bounds - see 
arrows. 
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Fig. 7 - J-i ntegral versus strain for various axi symmetric loadings 
applied to a circular crack. 
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- J-integral versus stress for various axi symmetric loadings 
applied to a circular crack. 
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Fig. 9 - Numerical results for a half-circular surface crack under 
uniaxial loading compared to two estimates. 
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Fig. 10 - Numerical results for a circumferential crack at a spherical 
void under uniaxial loading compared to estimates based on a 
circular crack. 
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Fig. 12 - Effect of geometric constraint on the plastic zone correction 
term for uniaxial loading of central cracks in plane stress 
and plane strain and of circular cracks. 
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Fig. 13 - Effect of multiaxial loading on the plastic zone correction 
term for circular cracks under axi symmetric loading. 
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